import gaussSeidel as GS
import LUdecomp as LU
A1=[[16.,3.],[7.,-11.]]#[[16.0,7.0],[3.0,-11.0]]
A2=[[2.,5.],[3.,7.]]
A3=[[2.,3.,-1.],[-1.,2.,3.],[0.,1.,2.]]#[[2,-1,0],[3,2,1],[-1,3,2]]
B1=B2=[11.0,13.0]
B3=[-1.,9.,5.]

x1=x2=x3=0
print '\nDoing x1'
x1=LU.doLUsolve(A1,B1)
#x1=GS.doMatrixSolve(A1,B1)
#print '\nDoing x2'
#x2=GS.doMatrixSolve(A2,B2)
print '\nDoing x3'
#x3=LU.doLUsolve(A3,B3)
#x3=GS.doMatrixSolve(A3,B3)
print ''
print 'x1:',x1
print 'Correct x1:',[0.8122,-0.6650]
#print x2
print x3
